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Abstract We consider a toy model of interacting extrovert and introvert agents in¬ 
troduced earlier by Liu et al [Europhys. Lett. 100 (2012) 66007]. The number of 
extroverts, and introverts is N each. At each time step, we select an agent at ran¬ 
dom, and allow her to modify her state. If an extrovert is selected, she adds a link 
at random to an unconnected introvert. If an introvert is selected, she removes one 
of her links. The set of N 2 links evolves in time, and may be considered as a set of 
Ising spins on an N x N square-grid with single-spin-flip dynamics. This dynamics 
satisfies detailed balance condition, and the probability of different spin configura¬ 
tions in the steady state can be determined exactly. The effective hamiltonian has 
long-range multi-spin couplings that depend on the row and column sums of spins. 
If the relative bias of choosing an extrovert over an introvert is varied, this system 
undergoes a phase transition from a state with very few links to one in which most 
links are occupied. We show that the behavior of the system can be determined ex¬ 
actly in the limit of large N. The behavior of large fluctuations in the total numer of 
links near the phase transition is determined. We also discuss two variations, called 
egalitarian and elitist agents, when the agents preferentially add or delete links to 
their least/ most-connected neighbor. These shows interesting cooperative behavior. 
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1 Introduction 

In recent years, there has been a lot of interest in the study of networks, many differ¬ 
ent types of networks have been studied: transportation networks like the railways or 
airlines networks 0, chemical reaction networks in a cell 0, power grids etc. a. 
A good review of the general theory may be found in j4j|5). In studies of networks in 
social sciences, some examples are the citation network of scientific publications |[6l, 
small-world networks 0. An important question that has attracted a lot of attention 
is the time evolution of social networks, in which the number of friends a particular 
agent interacts with evolves in time. In this context, a very interesting toy model was 
introduced by Liu et al J8], called the extreme introvert extrovert model. Numerical 
studies of this model revealed a phase transition, where the fractional number of 
links present undergoes a jump from a value near 0 to a value near 1, as a parameter 
in the model is varied continuously ei main I- Recently, we have shown that this 
model can be solved exactly in the limit of large number of agents m . This article 
is a brief pedagogical account of these results. For details, the reader is referred to 
the original publication. 


2 Definition of the Model 

We consider a group consisting of N introvert and N extrovert agents. The agents 
can add or delete links connecting them to other agents. In so doing, in our model, 
they do not have to get permission from the agent to whom the link is being added or 
from whom the link is being removed. The introverts are assumed to be comfortable 
with only a few links, and extroverts like many. In a general model, the nember of 
links an introvert likes to have is kj , and an extrovert likes kg, with kj < kg- We 
will consider the extreme case where k; = 0, and kg = Thus, and introvert does 
not like any links, and will delete a link, given an opportunity. On the other hand, 
an extrovert will try to add a link whenever possible. This model has bee called the 
eXtreme Introvert-Extrovert model(XIE) (8j . 101. 

A configuration of the system at any time is specified completely by an N x N ma¬ 
trix A whose entries are Ajj, with A ,/ = 1 , if there is link between the /-th introvert, 
and the /-th extrovert, and A ;/ = 0 otherwise. The total number of configurations is 
2 N ~ . The model undergoes a discrete-time Markovian evolution defined by the fol¬ 
lowing rules: At each time step, select an agent at random, and allow her to change 
the number of links connecting her to other agents. Any introvert has a probability 
(l+z)# , being selected, and an extrovert has a probability /^,\ N of being selected. 
Then, z < 1 corresponds to a bias favoring introverts, and z > 1 favors extroverts. 
If an introvert is selected, and has at least one link to an extrovert, she deletes one 
of the links connecting her to other agents at random. If she has no links, she does 
nothing, and the configuration remains unchanged. If an extrovert is selected, she 
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will add a link to one of introverts not already linked to her. If she already has all 
links present, she does nothing. 


3 Steady State of the XIE model 

We may think of the A 2 binary variables A J ; as Ising variables placed on ail A x A 
square grid, and then the update mle corresponds to single-spin-flip dynamics of the 
model. Note that our rules have the symmetry of changing introverts to extroverts, 
and Aij L-> 1 — A ,j , z 1 /z. This corresponds to the Ising model having spin-reversal 
symmetry with z+A 1/z. 

In general, given some such set of update rules, it is very difficult to determine the 
probabilities of different configurations in the steady state exactly. The remarkable 
result about the XIE model is that in this case, the probabilities of transition between 
configurations ^ and c €' satisfy the detailed balance condition, and one can write 
down the probability of different configurations in the steady state exactly. For the 
configuration 'A, in which the ;-th introvert has degree /?,, and the j- th extrovert has 
degree qj, the steady-state probability Prob* ( < A) has a very pleasing form (8) 

=^ w nw)fi(^-^)« (i) 

where f> (z) is a normalization constant. 

We may define the negative logarithm of this probability as the ‘energy’ of the 
configuration A’, giving 

= - £ log/?,! - £ log(N - qj)] - log(z) (2) 

i= 1 j =1 i 

We see that the effective hamiltonian has long-range couplings, and the energy 
of a configuration depends only on the row- and column- sums of the square array 
A. Also, the energy function is non-extensive: the energy of the configuration with 
all links absent varies as —A 2 log A. This non-extensivity causes no real problems, 
as all probabilities are well-defined, and averages of observables in steady state are 
well-behaved. 

Monte Carlo simulations of the XIE model have shown that, for large A, the 
system seems to undergo a transition from a few-links phase for z < 1 to a phase in 
which almost all links present for z > 1. In the few-links phase, the average number 
of links per agent remains finite, of order 1, even as A tends to infinity, with fixed 
Z < 1. Conversely, in the link-rich phase for z > 1, the average number of links per 
agent is nearly A, and the difference of this number from A remains finite, as A is 
increased to infinity. 

The fact that energy depends only on the row- and column- sums, and thus only 
on 2A variables, instead of the A 2 variable A, ; - explicitly, suggests that some kind of 
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mean-field treatment may be exact for this problem. This turns out to be true, as we 
proceed to show, but the treatment needs some care, as the number of variables, and 
hence also their conjugate fields, tends to infinity, in the thermodynamic limit. If 
the energy of the system depended only on the total number of links in the system, 
a single mean-field variable conjugate to the net magnetization would have been 
sufficient. 


4 Asymptotically exact perturbation theory 

We consider the low-density phase (z < 1) first. The case z > 1 is equivalent to this 
by the Ising symmetry discussed above. In the low density phase, the typical degree 
qj of the j -th extrovert is much less than N. Then we have (N — qj) ! « N\N~ q i . This 
suggests that we write for all q > 0 

(N-q)\/N\ =N~ q F(q,N) (3) 

with 

F(,,A0 = n( 1 -^) (4) 

For q F(q,N) is nearly equal to 1. Then, since W qj = we can write the 
effective Hamiltonian for the random XIE model as 

■ffleff ~ + ^int ( 5 ) 

where 

^0 = — £[ ln (/>;!)+/>< In -N\n{N\) (6) 

and 

J&, = -£lnF($;,A0 (7) 

j 

If we ignore the effect of the the ’’perturbation term ” different introverts 

are independent, and one can sum over states of each introvert separately. This gives 

f2 0 = (AHflcoof (8) 

with 

coo =Y J z k F{k,N E ) (9) 

k 

For large N, F tends to 1, and we get 

COo = 1 +z + z 2 +z 3 + ... = , (10) 

1 -z 
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which gives 

lo gQ(z) =N log (TV!) + N log (—!—) + .... (11) 

1 — z 

In a systematic perturbation theory, we need to determine the behavior of the 
steady state of the system under J%. This is easily done. In particular, we can de¬ 
termine the degree distribution of introverts and extroverts. It is easily seen that for 
large N, the probability that an introvert has degree r has an exponential distribution 
: Prob (introvert has degree r) = (1 — z)z r . Here, the /?,! factor in the weight of a 
configuration makes the usually expected Poisson form into a simple exponential. 
However, the degree distribution of the j-th extrovert is a sum of N mutually in¬ 
dependent variables ,4 ;/ , hence it remains a Poisson distribution. Clearly the mean 
degree of extroverts is same as the mean degree of introverts, so the Poisson dstri- 
bution has mean , which determines it completely. 

To lowest order in (1 /N), \ogF(q,N) = —q{q — 1)/(21V). Thus, while the inter¬ 
action hamilonian has different values for different configurations, and thus not a 
trivial c-number term, it is a sum of N different weakly correlated terms, and its 
mean is f)'( 1), and fluctuations about the mean are smaller. It is easily seen that they 
are (?(N ~'/ 2 ), giving 


logf2(z) = N log(AM) + N log (—'—) + 4?(1), for z < 1- (12) 

1 — z 

For z > 1, similar analysis, or the introvert-extrovert flip symmtery can be used 
to deduce that 

logf2(z) =N \og(N\)+N 2 \ogz+N log (-—) + ^(1), forz> 1. (13) 

1 - 1/z 


This is a remarkable result. Clearly, we get asymptotically exact result for 
log £2 (z) up to the linear order in N using the hamiltonian Mq. The effect of is 
only a term of f)’( 1) in £2 (z). In particular, in the I arge-/V limit, the density of links 
is 0 for z < 1, and 1 for z > 1. 

We note that these results are consistent with a scaling ansatz 


^(z)] 1/jV 

IV! 


= N a f(eN b ) 


(14) 


where z = exp(—e),a = b = 1/2, and /(0) = a finite constant, and f(x) is continu¬ 
ous at x = 0. For large positive x, f(x) ~ 1 /x. For z > 1, a similar form is obtained by 
the inversion symmetry, but the scaling functions for e > 0 and e < 0 are different, 
reflecting the ‘first-order nature’ of the transition. 
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5 Variants of XIE with preferential attachment 

It is interesting to consider some variations on this general theme. We consider the 
case when the agent does not choose which link to add (or delete) at random, as 
done in the XIE model in section 2, but decides on the basis of knowledge of degree 
of the receiving node. We consider two variations. 

Egalitarian agents: Here extrovert agents realize that the introverts regard links 
as burden, and attempts to distribute this burden as evenly as possible, and would 
add a link to the least connected introvert. In case Similarly, an introvert would cut a 
link to the most connected extrovert, as this action would make the other extroverts 
more equal. 

Elitist agents: Here, we consider the opposite extreme. In this case, an extrovert 
prefers the most ‘sociable’ introvert, and adds a link to the most connected of the 
available introverts. Similarly, an introvert cuts a link to the least connected avail¬ 
able extrovert. 

These variations have a strong effect on the degree distribution in the steady state. 
Let us discuss egalitarians first. Then, at any time, the degree distribution of an agent 
will have only two possible values: k or (k + 1), for some possibly time-dependent 
k. In the low density regime of this egalitarian XIE model, there are only a small 
number of contacts. It is easy to see that in the Iarge-/V limit, in the steady state, we 
have k = 0, and fractional number of introverts with exactly 1 contact is z. For the 
degree distribution of extroverts, it is easy to see that degree distribution is Poisson, 
with a mean value that increases with z. For z > 1, the only possible values of degree 
of an introvert are N — 1 and N. 

The behavior of the degree distribution at the phase transition point z = 1 is par¬ 
ticularly interesting. Here the fractional number of links can vary from 0 to 1, and 
the degree of an agent vary from 0 to N. However, the agents, by their cooperative 
behavior ensure that the inequality in the soceity always remains low. This is shown 
in Fig. [I], where we plot the time-dependent degree of two introverts, and two ex¬ 
troverts. While the actual degree varies in an unpredictable way, the four curves fall 
on top of each other, and are not distinguishable in the plot. 

We now discuss the elitists case. To recall, here the extrovert agents prefer to 
link to one most sociable of the introverts ( most connected), and introverts, to keep 
inversion symmetry, are assumed to delete their link to the least connected extrovert. 
This generates an interesting instability: Say, start with all links absent. Then, an 
extrovert will add a link to some randomly chosen introvert. Then all later extroverts 
selected at subsequent times will choose to link to the same agent. Her degree will 
tend to become N. Then the extroverts will choose a new ’’star’ to link to. The degree 
distribution of introverts becomes very unequal. The question what happens in the 
steady state? 

Interestingly, it turns out that in this case, for all z, the degree distribution is 
rather wide. Let us try to understand this. Clearly, argument sketched above needs 
elaboration. When a ‘star’ introvert’s degree becomes comparable to N, often an up¬ 
dating extrovert would be already connected to this introvert, and the time between 
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Fig. 1 Time trace of the degrees of two introvert k (black, blue), and of two extroverts, q (red, 
green) in a critical egalitarian XIE with N = 100. The graphs fall on top of each other, and are not 
distinguishable. Here the unit of t is a sweep. Taken from 03- 


changes of degree increases if the degree is closer to N. Meanwhile, other new stars 
are already in the making. This produces a broad distribution of degrees. If we in¬ 
spect the incidence matrix A at different times during the evolution of the system, 
we do not see much pattern in it directly. However, if we look at the same matrix 
after permuting the rows and columns to matrix so that both introverts and extro¬ 
verts are arranged according to their degrees in ascending order, we see much more 
structure. It is found, and easy to prove a posteriori, that in the sorted matrix, all l’s 
come in a single block with no holes, and similarly with zeroes [ Fig. 2], There is 
a single staircase-shaped interface that separates l’s and zeros, and this fluctuates 
with time. Clearly, if we have this property at one time, it will be preserved by the 
dynamical rules. 

The total number of accessible matrices in steady state is < (“^)(AH) 2 , which 
is much less than the number of possible 2 N ~. Thus, the configurations are much 
simpler to describe, when the permutation symmetry between agents is factored out, 
and one works with the equivalence classes of configurations under the symmetry 
group. As time evolves, the interface describing the configuration will evolve (we 
have to reorder the agents according to their degree distribution at each time). One 
can write the evolution rules in terms of the interface model. It can be seen that the 
interface will evolve by flipping at corners: height can change at the maxima, or 
the minima, at a rate proportional to length of side on the right [Fig. [3}. But this 
model seems difficult to solve exactly. On general grounds, since it is an interface 
fluctuation model where height can increase or decrease with equal probability, it 
may be expected to be in the universality class of Edwards-Wilkinson model, where 
the dynamical exponent is 2, i.e. relaxation time for an interface of length L varies 
as Lr. For the elitists case, this corresponds to relaxation time being approximately 
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Fig. 2 The incidence matrix after reordering. The 1 ’s and zeroes are separated by a single staircase¬ 
like interface. Here green = 1, yellow =0. 


@{\) sweeps (where one sweep of the system is L? attempts which updates each 
link on the average once), for large N. 

For different values of the z, the slope of the interface changes, but qualitative be¬ 
havior of relaxation remains the same. Hence, we find that the elitists self-organize 
into a critical state, where the degree distribution of agents is wide, for all z. 



Fig. 3 Evolution rules for the interface: height can change at the maxima (pink diamond), or the 
minima(blue diamond), at a rate proportional to length of side on the right. 


6 Summary and Concluding Remarks 

In this article, we discuss a simple toy model of a dynamical networks, where the 
state of different links keeps changing with time, but the system has a non-trivial 
steady state. As a function of the bias parameter z, the system undergoes a phase 
transition from a state with few links to a state with most links occupied. We showed 
that one can develop a perturbation theory, which becomes asymptotically exact for 
large N , at first order of perturbation. The corresponding state has non-trivial many- 
body correlations. We also discussed variations of the basic XIE model where the 
agents attach to one of the most- or least- connected agents, and showed that the 
steady state can have a wide degree distribution of degrees of agents. An interesting 
open problem is a theory to calculate the scaling function f(x) in Eq.(14) exactly. It 
is hoped that future research will lead to a better understanding of this model. 
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